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Abstract. We develop a model for the reflection and transmission of plane waves by 
an isotropic layer sandwiched between two uniaxial crystals of arbitrary orientation. 
£h In the laboratory frame, reflection and transmission coefficients corresponding to 

' the principal polarization directions in each crystal are given explicitly in terms 

£H . of the c-axis and propagation directions. The solution is found by first deriving 

O | explicit expressions for reflection and transmission amplitude coefficients for waves 

propagating from an arbitrarily oriented uniaxial anisotropic material into an isotropic 
material. By combining these results with Lekner's (1991) earlier treatment of waves 
propagating from isotropic media to anisotropic media and employing a matrix method 
we determine a solution to the general form of the multiple reflection case. The example 
I/"") \ system of a wetted interface between two ice crystals is used to contextualize the results. 
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1. Introduction 



Since the time of Descartes scientists have taken an interest in the physical properties of 
ice (Dash et al. 2006). For example, Tyndall (1856) observed harvested ice in some detail, 
but it was not for almost another hundred years, when Nakaya (1954) began cataloging 
observations of snowflakes, that well-controlled laboratory growth experiments on 
single crystals began. More recently, motivated by the fact that commonly occurring 
environmental temperatures span the triple point of ice, a strong understanding of the 
thermodynamics and phase behavior of polycrystalline ice near its melting temperature 
has been developed. It is known that near the melting temperature, an interconnected 
network of liquid water exists within the polycrystalline solid (Nye & Frank 1973). Under 
hydrostatic conditions, veins of water separate the boundaries between three crystals 
and join in nodes where four grains meet. The Gibbs-Thomson and impurity effects are 
responsible for the presence of this liquid (Dash et al. 2006), which is observed using 
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optical microscopy techniques (Mader 1992, Walford et al. 1987). Additional water 
structures, such as water lenses, are observed in the presence of non-hydrostatic stresses 
(Nye 1991). Less well understood is what happens away from these junctions in the 
planar interface between two single crystals. While disorder is expected at the molecular 
scale, it is predicted that a dopant, such as salt, can induce the formation of a thick (> 10 
nm), essentially bulk, water film (Benatov & Wettlaufer 2004). This phenomenon of 
interfacial melting could have important implications for ice's electrical and mechanical 
properties and impurity redistribution in glaciers and polar ice (e.g. Rempel et al. 2002). 
More generally, it may occur in other polycrystalline materials. Recent predictions 
of interfacial melting at ice grain boundaries have motivated an experimental search 
(Thomson et al. 2005) to detect the water layer using an optical reflection technique. 
This has led us to consider the theoretical formalism for wave reflection and transmission 
in an anisotropic/isotropic/anisotropic layered system; specifically when the anisotropic 
media are uniaxial crystals. 

Uniaxial crystals are scientifically well studied materials owing to their ubiquity 
in nature and their many technical applications, including use as elements in optical 
systems. Theoretical treatments of light propagation in these crystals have focused 
on reflection from surfaces, and the internal propagation through layered structures. 
Previous studies have used 4 x 4 or 2 x 2 matrix methods to solve the general problem of 
light propagation through birefringent networks, where solutions are given in principal 
axes coordinate frames (Yeh 1979, Yeh 1982). These treatments leave the reader to 
solve the eigenvalue problem associated with transforming to a laboratory coordinate 
frame of experimental relevance. In other studies only special c-axis orientations are 
considered (e.g. Yeh 1982) or multiple internal reflections are ignored (e.g. Gu & 
Yeh 1993). Still other methods find solutions at a single interface and are, due to their 
form, difficult to extrapolate to multiple interfaces or multiple reflections (Stamnes & 
Sherman 1977, Zhang & Caulfield 1996). These existing theoretical studies are not 
ideally suited to experimental applications; to be of utility solutions must be valid for 
arbitrary crystallographic orientations, incidence angle, and propagation direction, all 
measured in the laboratory frame. 

Our study approaches the problem from the perspective of the specific experimental 
setting described above; a wetted interface between two uniaxial crystals of arbitrary 
orientation. We begin by revisiting the problem of plane wave propagation in a uniaxial 
crystal using modal decomposition, the approach by which Lekner (1991) determined 
the reflection and transmission amplitudes for a plane wave entering an anisotropic 
medium from an isotropic medium. The relevance of Lekner's (1991) important earlier 
work to our study requires that we begin by reviewing his results in some detail. Here we 
generalize those results in order to analyze the reverse situation; where the wave passes 
from an anisotropic to an isotropic region. We explicitly determine the reflection and 
refraction amplitude coefficients in terms of the orientation of the c-axis with respect to 
the laboratory axes and the optical constants of the materials. Consequently, we are able 
to solve for all of the relevant amplitude coefficients associated with an isotropic layer 
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sandwiched between uniaxial crystals. This enables us to construct a matrix method to 
model the Fabry- Perot effect of multiple reflections from the isotropic layer. Extensions 
of the theory may also be applicable to more general birefringent systems, but here it 
is of particular interest at the interface between two grains in water ice. To illustrate 
this we present clear examples of how the generalized theory can be used in comparison 
with light reflection experiments. 

2. Anisotropic Optical Theory 

The problem of interest, illustrated schematically in figure [U(a), is that of plane wave 
propagation in a three layer system, an isotropic layer (Z 2 ), bounded by uniaxial crystals 
{li and Z 3 ). Within the isotropic layer a wave's polarization can be decomposed in 
the usual way to be parallel (p) and perpendicular (s) to the plane of incidence. In 
the anisotropic material the principal components are parallel and perpendicular to the 
optical axis of the material; these are referred to as the extraordinary (e) and ordinary (o) 
modes. To characterize the system completely we describe the plane wave propagation 
within each layer, in addition to the reflection and refraction at the boundaries between 
the media. 

In the laboratory frame of reference the reflecting surfaces are xy planes and z is the 
normal; the zx plane is chosen as the plane of incidence. The electric field is denoted 
E = [E x , E y , E z ]e^ qz+Kx ~ ujt \ No y dependence exists due to translational symmetry 
in the y direction. Continuity of the tangential component of E demands that K is 
common in all media while the normal component of the wave vector (q) will depend 
on the state of polarization, propagation direction and specific medium of propagation. 
For example the value of q corresponding to an ordinary ray propagating in the -z 
direction is denoted by q~ . The angular frequency is lu, thereby defining a wavevector 
k = uj/c. Referring to figure d^b) for waves incident from an anisotropic media onto an 
anisotropic/isotropic boundary reflected o and e waves and colinear, transmitted s and 
p waves result. Conversely (Figure [D(c)), for waves incident from an isotropic media onto 
an isotropic/anisotropic boundary o and e waves are transmitted and colinear s and p 
waves are reflected. For the three layer system, a given K determines unique qi SO and 
8 iso values within the sandwiched isotropic layer. Alternatively, if one was interested in 
a single angle of incidence in the uniaxial crystal the problem could be solved iteratively 
using different K values for the o and e polarizations. 

2.1. Plane Wave Propagation within a Uniaxial Crystal 

Lekner (1991) employed a normal mode analysis to study wave propagation in uniaxial 
crystals. Using an orthogonal coordinate transformation, expressed in terms of direction 
cosines, he first expressed the dielectric tensor of a uniaxial crystal in the laboratory 
frame of reference. He then explicitly determined the ordinary (E°) and extraordinary 
(E e ) electric field vectors, in addition to the wave vector's z-components for propagation 
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Figure 1. (a) The general geometry of the three layer system. Plane wave propagation 
is illustrated with rays, enlarged schematics of the boxed interface regions are shown 
in (b) and (c) and in 1% an example crystal is inset. The c-axis direction of the example 
crystal is labeled with respect to the laboratory frame of reference. Here 9 a is the angle 
between c and x, Op is the angle between c and y, and 1 is the angle between c and 
z. (b) Schematic of an wave incident on an anisotropic/isotropic boundary, where 
K is the preserved tangential component of all wave vectors, g's are the normal (z) 
components of the wave vectors, and 8i SO is the angle of the transmitted wave vector. 
An analogous schematic of an incident e wave could be drawn, (c) An s or p wave, 
they are co-linear, incident on the isotropic/ anisotropic interface. 



within the crystal, in terms of K, k, and the c-axis orientation, specified by the unit 
vector c = [a,/3, 7]; where a = cos 6^, (3 = cos 6^, 7 = cos# 7 (see figure [Qa)). 

E 0± = Nf[-Pqi, aqf - jK, f3K], (1) 
E e± = NtW D - wtK, (3e k 2 , 7M 2 - qf) ~ aqfK). (2) 

Here q Q and q e are the normal components of the ordinary and extraordinary wave 
vectors, e Q = n 2 Q is the ordinary dielectric constant, and N Q and N e are normalization 
constants. The expressions for the z-component of the wave vectors are 

q± = ±^e D k 2 - K 2 and (3) 
, _ ±VD-a 1 KAe 

9e e + 1 2 Ae ' [) 

where in all cases the signs (±) correspond to the direction of beam propagation with 
respect to the z-axis. The quantity D, in qf, is given by 

D = e [e e (e + 7 2 Ae)fc 2 - (e e - 1 Ae)K\ (5) 

where e e is the extraodinary dielectric constant and Ae = e e — e . 

Lekner's (1991) analysis is general to plane wave propagation within uniaxial 
materials and provides a foundation for investigating reflection and refraction at 
interfaces with such materials. 
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2.2. Lekner's Amplitude Coefficients 

Lekner (1991) goes on to calculate reflection and transmission coefficients for s and p 
waves incident from an isotropic media onto an isotropic/anistropic interface. Because 
the tangential components are preserved across the boundary he focuses on the z- 
dependence of the electric field. The z-dependencies of the incident, reflected, and 
transmitted electric fields for the s polarization are, 

incident: E in = e^-*[0, 1, 0], 

reflected: E ref = r sp e iq ™° z [cos 6 iso , 0, - sin 6 iso ] + r ss e^° z [0, 1, 0], (6) 
transmitted: E tr = t so e iq ° z [E°\ E° + , E° + ] + t se e iqtz [E e x + , E e y + , E e z + }. 

At the interface the electromagnetic waves are subject to the following boundary 
conditions implied by Maxwell's equations; continuity of E x , E y , dE x /dz — iKE z , 
and dEy/dz, where the subscripts x,y,z refer to the vector components. Applying 
these boundary conditions at the reflecting plane (z = 0) leads to four equations that 
Lekner (1991) solved for the four unknown intensity coefficients: r ss ,r sp ,t so , and t se . 
In Appendix A ( 1A.ll) -( lA~4l) we summarize those results with small changes correcting 
apparent typographical errors in the original publication. The incident, reflected, and 
transmitted p waves are 

incident: E m = e iq ™° z [cos 9 iso , 0, - sin 9 iso ] , 

reflected: E ref = r pp e iq ™° z [cos6 iso , 0, - sm6 iso ] + 7>e^°*[0, 1, 0], (7) 
transmitted: E tr = t po e iq ° z [E°\ E° + , E° z "] + t pe e iqt z [E e x \ E° + , E e z + }. 

Again, application of the four boundary conditions at the interface leads to four 
amplitude coefficients, r pp ,r ps ,t po , and t pe , also presented in Appendix A flA.8j) - flA.llj) . 

When considering these results for transmission from isotropic media to anisotropic 
media it is critical to realize that within the anisotropic material the ray direction (i.e., 
the Poynting vector) for the extraordinary mode differs from the wave vector direction. 
This subtlety must be recognized to verify simple test cases of reflection and refraction. 
For more discussion regarding ray direction see Lekner (1991). 

2.3. Anisotropic to Isotropic Interface 

We analyze the reverse incidence, when o and e waves are incident from an anisotropic 
media onto an anisotropic/isotropic interface, in a similar manner. For the ordinary 
wave the z-dependencies of the electric fields are 

incident: E m = e iq ° z [E° + , E° + , E° z + ] , 

reflected: E ref = r 00 e iq ° z [E° x ~,Ef,Ef] + r oe e iq ° z [Ef,Ef,Ef], (8) 

transmitted: E tr = t os e iq ™° z [0, 1, 0] + t op e iq ™° z [cos6 is0} 0, - sm6 iso \. 

Applying the boundary conditions to the ordinary wave at the reflecting plane (z = 0) 
the above general expressions for the incident, reflected, and transmitted waves (IE]) yield 



Light scattering from an isotropic layer between uniaxial crystals 



6 



a system of four equations; 

E° + + r 00 E°~ + r oe El - t op cos 6 iso = 0, (9) 

E° y + + r 00 Ef + r oe Ef - t os = 0, (10) 
(lo E °x + + % r ooE°~ + q~r oe E e ~ - qf so t op cos6 iso - 

K(E° Z + + r 00 E°~ + r oe Ef + t op sin 6 iso ) = 0, (11) 

q+Ef + q-r 00 E° y + q~ e r oe E e y - q+ so t os = 0, (12) 

and four unknown amplitude coefficients r QO ,r oe , t os and t op for the ordinary wave. 
Solving this homogeneous system of equations for the unknown amplitude coefficients 
provides expressions for r OQ ,r oe , t os and t op shown in their complete form in Appendix B 
(JEQ-flB3]). 

The extraordinary wave's intensity coefficients can be found in a manner analogous 
to those for the ordinary wave. Again we begin with expressions for the z-dependence 
of the extraordinary electric field: 

incident: E in = e iqtz [Ef , E e y + , E?] , 

reflected: E ref = r ee e^ z [E e x ~ , Ef , Ef } + r eo e^ z [Ef , Ef, Ef ],(13) 

transmitted: E tr = t es e ig ™° z [0, 1, 0] + t ep e ig ™° z [cos 6 iso , 0, - sm6 iso \. 

We find solutions for the intensity coefficients of an incident extraordinary beam (r ee ,r eo , 
t es and t ep ) as we did previously; see ( 1B-5I) - (1B-8I) . Thus, the magnitudes of the derived 



amplitude coefficients ( 1B-1I) - (1B-8I) are fully determined by completing the normal mode 
analysis substitutions ([I])-©. 

A limiting case provides some verification of the now explicit amplitude coefficients 
for an incident o wave (1B-1|) -( |B^4|) presented in Appendix B. If the c-axis of the crystal 
is in the plane of incidence and the incident o wave is entirely perpendicularly polarized 
with respect to the plane of incidence (i.e. [E°,E°,E°] = [0,1,0]) the reflection and 
transmission amplitude coefficients reduce to the Fresnel equations for perpendicular 
polarization (e.g. Born & Wolf 1965): 



qj - qJso Qo - Qls 



Visa - % <itso + Qt 



and (14) 



q%so Qo Qiso ' Qo 

remembering g+ = — q~ . As expected the other coefficients (r oe , t op ) vanish. In contrast, 
the extraordinary amplitude coefficients can be modeled using isotropic theory for the p 
polarization as long as it is recognized that this introduces an effective index of refraction 
that is a function of incident angle (Born & Wolf 1965), 

n e ff = =. (16) 

y/ n"l sin Q\ + n\ cos 9[ 

Here Q\ is the wavevector incident angle on the boundary and n Q and n e are the ordinary 
and extraordinary indices of refraction. It also must be noted that within the crystal 
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Figure 2. Extraordinary reflection and transmission coefficients for a beam incident 
on the basal plane ([9 a , 9p, 6~\ = [90°, 90°, 180°]) compared with isotropic theory. Here 
the Fresnel equations for p polarization solutions are used with n e f f (|16|) . Circles are 
isotropic theory and the points (inside circles) are the full theory. The line at zero 
represents the cross term coefficients, which for this crystal orientation are always 
zero. The physical constants used are n a = 1.1, n e = 1.2 and rii SO = 1.33. 



the wave vector and electric field are not necessarily perpendicular. Therefore, the angle 
used to compute the Fresnel coefficients must be that of the Poynting vector (the ray 
direction), while Snell's law must be solved using the wave vector direction. Figure [2] 
illustrates the agreement between the coefficients for an extraordinary beam incident on 
a basal plane and the isotropic Fresnel equations. 

Now that we have investigated plane wave propagation within each region and 
across each boundary, individually; we return to the three layer system. By collecting 
the expressions for each of the relevant reflection and transmission amplitude coefficients 
it is possible to construct a matrix formulation for the propagation of light through the 
uniaxial network. 

3. Matrix Method 

Similar to the Jones matrix formalism, 2x2 matrices can be used to describe reflection 
and refraction at interfaces with uniaxial materials (e.g. Yeh 1982, Abdulhalim 1999). 
Rather than rotation matrices, as in the Jones formulation, here the matrix elements 
are the relevant amplitude coefficients. The diagonal elements represent reflection or 
refraction of like polarization, while the off-diagonal elements represent the mixing 
of polarization states. Each interface is represented by two independent matrices, 
one representing reflection, the other refraction. For example, reflection at an 
anisotropic/isotropic interface can be written in matrix form as 

flA (17) 
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where both the reflected wave, E lr = (E° r , Ef r ), and the incident wave, E; = (E°, Ef), 
will have ordinary and extraordinary components. The reflection matrix Ri is composed 
of the amplitude coefficients representing the interface which incorporate the properties 
of the anisotropic material. The phase shift acquired by waves that travel some distance 
through a uniaxial material can be accounted for using diagonal propagation matrices, 
such as 

A)- 

The phase factors for the o and e waves are given by 5 a = A oA / 1 q 2 + K 2 and 5 e = 
K e ^J q 2 + K 2 , where A Q and A e are their respective path lengths within the crystal. 

This matrix approach for treating different interfaces substantially simplifies the 
analyses of light scattering associated with layered materials. However, it is important 
to note, when modeling multiple layers and/or reflections careful attention must be paid 
to maintaining a consistent coordinate system. In the next section we use this approach 
to examine the specific example of an isotropic layer sandwiched between anisotropic 
layers. 



4. Anisotropic/Isotropic/ Anisotropic Layering - Application to Grain 
Boundaries 

As discussed in the introduction, a primary motivation of this study is to better 
understand the phase behavior of ice and other polycrystalline materials near their 
melting temperatures. One proposed experimental method for characterizing the 
grain boundary is to measure a reflected laser beam's intensity as a function of the 
thermodynamic variables: temperature, crystal orientation, and impurity concentration 
(Thomson et al. 2005). However, data gathered from such an experiment can only 
be interpreted accurately with a theoretical model that includes the anisotropy of the 
system. The combination of our results with those of Lekner (1991) leads to precisely 
that type of model; of an isotropic layer sandwiched between uniaxial crystals (Figure[3]). 
While here we focus on ice and water to make a connection to our experiment, the 
theoretical model applies to any such geometry and anisotropy. 

First we ignore any optical path length within the ice crystals, and simply consider 
the problem as if the incident beam were generated and the reflected beam observed 
at the interface of l\ and I2 (Figure [3]). Propagation away from the interface can be 
accounted for using (fl8|) . Following the matrix approach of Section [3] the transmission 
and reflection coefficients for each interface can be formulated; 




Light scattering from an isotropic layer between uniaxial crystals 



9 




Figure 3. Schematic of multiple reflections off of a grain boundary interface. Different 
interfaces are labeled with the relevant amplitude coefficient matrices for reflection and 
transmission. The isotropic film (1 2) between the crystals, l\ and Z3, has thickness d. 



and R\ remains as previously denned by (TT7l) . Primed amplitude coefficients are used 
to denote those associated with the h/h boundary. For now we ignore the transmitted 
signal, which does not apply to our experiment. For this model we define the page to be 
in the xz plane with positive up (z) and to the right (x), and are careful to insure that 
the polarization vector direction is consistent in each layer. Therefore, in the unprimed 
amplitude coefficient solutions of T 3 and _R 3 , obtained by solving ( 1A.ll) -( lA7lll . 9 iso is 



substituted with —9 iso , and the negative normal mode solutions from CQ)-(JI]) are used. 

The waves of experimental interest are the initial and each subsequent reflection 
(Ei r , E 2r , E 3r ...etc), and can be expressed as a function of the original incoming field 
as was done in ffl7|) for the E lr term. Writing down the first few reflections the pattern 
is evident; 

E 2r = T 3 R 2 T 1 E i e i{wt - s) , (21) 
E 3r = T 3 i? 2 i? 3J R 2 T 1 E,e^- 2<5 \ (22) 
E 4r = T 3 i? 2 i? 3 i? 2 i? 3 i? 2 T 1 E i e^"' i 3S \ (23) 

These terms include the additional phase contributions arising from the optical path 
length of each reflection internal to the grain boundary, 5 = 2 kni SO d cos 6i SO (Figure [3]). 
Because the multiple reflections occur within the isotropic medium, the path length 
depends upon only one angle. Anisotropy within the intervening layer would further 
complicate the situation by introducing a second angle of transmission resulting in 
multiple possible paths within the layer. In the limit of the superposition of a large 
number of such reflections, the total reflected field (E*°*) becomes 

E*°* = E lr + E 2r + E 3r + + E nr (24) 



= {R x + T 3 R 2 R n - 1 e- in5 T 1 )E i e iuJt (25) 

n=l 

where, R = -R 3 -R 2 . This expression contains a geometric series of matrices (e.g. 
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Strang 1993) and as such can be rewritten, using the identity matrix /, as 

E*°* = [#! + T 3 R 2 e- l5 (I - Re- i5 Y 1 T 1 ]E i e iwt . (26) 

This substitution is valid as long as the absolute values of the eigenvalues of R are less 
than one. The limiting case of total reflection with no anisotropy, R = (o i)> illustrates 
that any other situation will lead to smaller eigenvalues; less than one. The bracketed 
expression in (1261) will remain valid for reflection in any geometry that includes an 
isotropic sandwich. We call this the n-reflection matrix (M^ e ^), 

M ref _ ( [ r o e t2S + {Utpo + Cttso ~ r)ir 00 )e lS + fyfatop + °i] 3 r 00 + V^tos)] 
\ [r e oe m + {(itpo + (, 2 t so ~ Vl r eo)e i5 + rj2(Gtep + m r eo + V^es)] 

[r oe e l2S + (( 4 t P e + Cbtse - r] 1 r oe )e lS + r] 2 (( 6 t op + n 3 r oe + ( 7 t os )} \ 1 

[r e ee m + (Cltpe + Cltse ~ V^ee)^ 5 + m{(4e P + ^ ee + C,-jt es )\ J (e l2S - ^ + rj 2 T] 3 ) ' 

(27) 

In (T27j) the variable quantities within the matrix elements are given in Appendix C. In 
an experimental system with a beam propagating through a crystal, the propagation 
matrices (ITS1) associated with the optical path length will modify the final result. 

This is the extent to which we can easily pursue the problem algebraically. 
Calculating an experimentally useful value, such as the reflected flux density (i.e. I r = 
E*°* • E*°'*/2), involves taking the complex conjugate of E*°* and solving for the reflected 
intensity by working through the full algebraic expressions. This quickly becomes quite 
cumbersome as can be seen from the full expression in (!27j) for M™* but is easily 
facilitated by a symbolic or matrix based mathematical computer interface. Examples 
of intensity ratio (i.e. Ir = E* ot ■ E*°**/(Ej • E*)) are plotted for reflection from a basal 
plane, as a function of incidence angle (Figure HJ) and as a function of grain boundary 
film thickness (Figure [5]). For other orientations plots of intensity ratio illustrate some 
interesting characteristics. For special crystallographic orientations (Figures EEf6]) no 
polarization mixing takes place. However, when a particular polarization state within 
the crystal is parallel with the p polarization in the isotropic layer, a Brewster-like 
angle exists for the system. In figure H] this Brewster angle is present for the e polarized 
wave, and conversely, in figure [6] the o polarized incident beam has an angle of zero 
reflection. For systems with less crystallographic symmetry (Figure [7]) the mixing of 
polarization states is clearly important. An examination of intensity ratio as a function 
of c-axis orientation (Figure [8]) illustrates the relative importance of polarization mixing 
and points where symmetries preclude coupling. Theoretical curves for experimentally 
measured crystallography (Figure [9]) show substantial polarization mixing, yet for 
certain orientations and polarizations effective Brewster-like angles emerge. 

A nearly identical analysis can be done in order to calculate the theoretically 
transmitted fields and intensity ratios associated with an isotropic layer sandwiched 
between anisotropic media. The infinite sum is simply rewritten in terms of the 
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Figure 4. Theory for reflected intensity ratio (Ir) from the basal planes, [6 a , 6p, 6> 7 ] = 
[90°,90 o ,0°] of each crystal, in an ice-water-ice sandwich. Intensity ratio is plotted 
for varying water layer thicknesses, d = 1, 10, 100 nm. In both cases solid curves 
correspond with reflected o polarization and dashed curves with the e polarization. 
Curves asymptote at the approximate critical angle for total reflection within the 
water layer (« 80°). (a) The incident beam is o polarized [E°,Ef] = [1,0]. Note that 
in this geometry, for purely ordinary incidence, the extraordinary reflected component 
is always zero, (b) The incident beam is e polarized [Ef,Ef] = [0, 1]. Again there is 
no polarization mixing. 




100 200 300 400 500 600 

d(nm) 



Figure 5. Theory for reflected intensity ratio (In) from the basal planes, [9 a , Op, 6* 7 ] = 
[90,90,0] of each crystal, in an ice-water-ice sandwich as a function of water layer 
thickness d, at an incidence angle of 55°. The solid curve corresponds to the reflection 
of an incident wave of purely o polarization [E°,Ef] = [1,0] and the dashed curve to 
the reflection of an incident wave of e polarization [E°, Ef] = [0, 1]. 
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Figure 6. Theory for reflected intensity ratio (Ir) from an ice- water-ice sandwich. 
[6£, 6^,6>i] = [90°,0°,90°] and [0£,0|,0^ = [90°, 90°, 0°]. Intensity ratio is plotted 
for varying water layer thicknesses, d — 1,10,100 ran. In both cases solid curves 
correspond with reflected o polarization and dashed curves with the e polarization, 
(a) The incident beam is o polarized [Ef,Ef] = [1,0]. (b) The incident beam is e 
polarized [Ef, Ef] — [0, 1]. No polarization mixing occurs in either case. 




Figure 7. Theory for reflected intensity ratio (In) from an ice- water-ice sandwich. 
[0*,0j3,0}] = [45°, 45°, 90°] and [0l,6j,6*] = [90°, 90°, 0°]. Intensity ratio is plotted 
for varying water layer thicknesses, d = 1,100,200 nm. In both cases solid curves 
correspond with reflected o polarization and dashed curves with the e polarization, 
(a) The incident beam is o polarized [E°,Ef] = [1,0]. (b) The incident beam is e 
polarized [Ef, Ef] = [0, 1]. Polarization mixing occurs in both cases. 
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(a) (b) 

Figure 8. Theory for reflected intensity ratio (Ir) from an ice-water-ice sandwich for 
multiple values of 0i SO . The c-axis orientation of l\ is rotated about the z-axis. The 
angle, <pg, is the projection of c onto the xy plane with 6* = 45° held constant. Water 
layer thickness is also a constant, d = 100 nm, and [8^,8^,8^] = [90°, 90°, 0°]. Again 
solid curves correspond with reflected o polarization and dashed curves with the e 
polarization, (a) The incident beam is o polarized [E°,Ef] = [1,0]. (b) The incident 
beam is e polarized [E°,Ef] = [0,1]. In both cases the curves for <pg = 180° — 360° 
are symmetric with what is shown. 




Figure 9. Theory for reflected intensity ratio (In) from an ice-water-ice sandwich 
grown in an experimental ice growth cell; [8^,81,8^] = [103°, 131°, 45°] and 
[0l,9p,9j] = [89.7°, 88.5°, 1.5°]. Intensity ratio is plotted for varying water layer 
thicknesses, d = 1, 100,200 nm. In both cases solid curves correspond with reflected 
o polarization and dashed curves with the e polarization, (a) The incident beam is o 
polarized [E?,Ef] = [1,0]. (b) The incident beam is e polarized [E°,Ef] = [0,1]. In 
both cases the mixing of polarization states is clear. 
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transmitted waves, 

oo 

E Tt = (T 2 Y,(Re~ lS ) n Ti)E t e^ (28) 

n=0 

= [T 2 (I - ^e- w )- 1 Ti]E i e <wt , (29) 
where we can label the bracketed portion of the expression the n-transmission matrix, 



M tr = ( [to P t' po + t os t>Je m + [ V5 t' po + ifctfje* 
\t Kn tL + t e J'Je i2S + \mt'„ 



[t ep t' po + t es t'Je m + [ V7 t' po + m t'Je iS \ 1 
[tepfpe + test' se ]e m + W pe + jfefje" j (e* 25 - Vl e^ + mm ) ' 

Again the variable quantities witidhin the matrix elements are given in Appendix C. 
Ensuing calculations follow in analogy to what we have previously shown for reflection. 
Whereas for our experimental system this is not of interest, there may be experimental 
systems in which the transmitted wave would be the appropriate observable. 



5. Conclusions 



We have derived explicit expressions for the reflection and refraction amplitude 
coefficients for ordinary and extraordinary polarized electromagnetic waves incident 
upon an interface between a uniaxially anisotropic and an isotropic material. The 
orientation of the optical axes may be arbitrary with respect to a laboratory frame of 
reference. Furthermore, the formulae are valid for the full range of incidence angle. 

Combining these results with Lekner's (1991) earlier work allows us to model a three 
layer system of arbitrarily oriented uniaxial crystals sandwiching an isotropic layer. A 
Jones-like matrix formulation is used to treat the electromagnetic wave propagation 
and the reflection and transmission matrices for an isotropic layer sandwiched between 
anisotropic materials are explicitly determined. Light scattering from the interfaces 
between two grains in polycrystalline ice is an area with wide ranging implications 
in astrophysical and geophysical settings, but also serves as an ideal, transparent 
analogue for many materials (Dash et al. 2006). Comparable systems may be present 
in layered ceramics (Luo & Chiang 2008), biological structures (Parsegian 2006), or 
experimental tests of the theory of dispersion forces in such layered geometries (van 
Benthem et al. 2006). Presently we are engaged in a long term experimental test of 
the theoretical framework described here (Thomson et al. 2005). It is hoped that, in 
consequence, those interested in the structure of grain boundaries in other systems can 
make use of both the framework laid out in this paper and the incipient experimental 
findings. 
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Appendix A 

Here we summarize Lekner's (1991) results for the reflection and transmission coefficients 
for s and p waves incident onto an isotropic/anisotropic interface. We have made small 
changes where apparent typographical errors were made in the original manuscript. 

_ (q+ - q^AEf - (q+ so - q+)B t E° + 



D 



i 

r. 



(A.1) 



r _ 2(g£ o cos0 iao + K^A^jAiEf - BjEl 
r sp — p. l A - z J 

tso = ~ 2q ™° Bl and (A.3) 

tse = (A.4) 

where 0i SO is the incident angle in the isotropic layer. The subscripts p,s,o and e denote 
the particular amplitude coefficient, the q's are the wave vector's z components, and 

M = (qt + qtso + K tan 6 lso )Ef - KE° Z + , (A.5) 
Bi = {qt + qfso + K tan 9 lso )Ef - KE e z + , (A.6) 

A = (qtso + ?M< - ItiL + ( A -7) 

Equation ()A.4j) has been altered after Lekner (1992) where he pointed out a misprint of 



the sign in Lekner (1991). The coefficients for the incident p wave are: 

2 (qtso cos °iso + K sin °iso)Fl < / a 

r pp — F> a ' 

Di cos iso 

2(g+ Q cos 6 iso + K sin 9 iso ) (qj - q+)E° + E? 
r ps = n — , (A.9) 

2(qtso cos Oiso + K sin e iso )(qt so + q+)E° + 
t po = u — and (A. 10) 

-2(qt so cosd lso + Ksme tso )(qt so + qt)E° + 

t pe = =: — , (A.ll) 



where F L = [(g+ D + q+)E° + Ef - (q+ so + qt)E° + E e x + }. Here (jA8]) differs from Lekner 
(1991), where we infer there was a typographical error. 

Appendix B 

The full expressions for the amplitude coefficients associated with the ordinary wave at 
an interface with an isotropic material are, 



EfjqL - qt)[Ef(A - q-) + KEf] - Efj^ - q~)[E° + (A - g+) + KEf] 

~ g , (B-l) 
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where, B = Ef (qf - qf [Ef (A -qf + KEf] - Ef (q+, - qf[Ef (A -qf + KEf] 
wdA = q++Ktfm6+ 



E°i(qL ~ qf[E°f(A - qf + KEf] - Ef (qf - qf[Ef (A - qf + KEf] 
B • (B_2) 



tos 

B- 1 [Ef {Ef (g+ - qf [Ef (A - qf + KEf] - Ef (g+ - qf [Ef (A - qf + KEf]}+ 

EfEf(q; - q-)[Ef{A - qf) + KEf]] (B-3) 

top 

(B cos e^f^Efiq^ - q-)[-EfEf{qt - qf + K(Ef E°f - EfEf)] + 
Efiqto ~ q- e )[EfEf(qt - qf + K(Ef E°f - Ef E°f)] + 

Ef(qL ~ qt)[EfEf{q- - qf + K(Ef Ef - Ef Ef)]} (B-4) 

For an incident wave polarized in the extraordinary direction the coefficients have 
a similar form; 



EfjqL - qo)[Ef(A - qf + KEf] - Ef(qf - qf)[Ef (A - q~) + KEf] 

— ^— 2 (B-5) 



f'eo 

Ef (qtso ~ qf) [Ef (A - qf + KEf] - Ef (qf - qf [Ef (A - qf + KEf] 

B 



(B-6) 



B- 1 [Ef {Ef (g+ - qf[Ef (A - qf + KEf] - Ef (g+ - qf [Ef (A - qf + KEf]}+ 

EfEf{q- - q-)[Ef{A - qf) + KEf]] (B-7) 



{Bcos6 lso )-fEf{qi so - qf[EfEf{qt - qf + K{Ef Ef - EfEf)] + 
Ef{qto ~ qo)[-EfEf{qt - qf + K(Ef Ef - Ef Ef)] + 

Ef(qL ~ qf[EfEf{q~ - qf + K(Ef Ef - Ef Ef)]} (B-8) 



Light scattering from an isotropic layer between uniaxial crystals 



17 



Appendix C 

Although any computation utilizing the presented theory is most efficiently done 
utilizing linear algebra, the variable quantities that compose the elements of the n- 
reflection matrices {M™ and Mff) are presented here for completeness. 



Vi = Rn + R22 


(c-r 


11 11 

rir. — r r — TV 
1^ ps sp pp ss 


(C-2) 


Tin = TV — V T 
'IS — ' ps 1 sp 1 pp 1 ss 


(C-3; 


= Tpptso ^sptpo 


(C-4; 


V5 = Rutos ~ R22top 


(C-5) 


V& = R2\top — Rlltos 


(c-6; 


7] 7 = Rutes — R22t e p 


(C-7) 


7]s = R2\t ep — Rlltes 


(c-8; 


Cl = r pptep + r S ptes 


(C-9; 


C2 = r ps tep + r sstes 


(C-10; 


C3 = T sstpo ^pstso 


(c-ir 


= r pptop + r S ptos 


(c-12; 


C5 = r ps top + r sstos 


(c-13; 


C(> = T sstpe ^pstse 


(c-w; 


C7 = ^pptse ^sptpe- 


(c-15; 



Where the subscripted i?'s are the matrix elements of the previously defined R, 

(ryt ,y I <y> ^ <y <y> ry^ I ry <y ^ 

rppfpp + r r sp Tpp r + r r ss . 

ry' ry I ry' ry ry ry' I ry' ry ' 

1 pp 1 ps T ' ps i ss 1 ps 1 sp ~T 1 ss i s 
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